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ABSTRACT 


A reliable system for performing a Fourier analysis of a 
reactor noise signal has been constructed and utilized. Both 
auto- and cross-correlation techniques were used to study. 
noise signals from the MITR hohlraum and the exponential 
lattice. As predicted by theory, the efficiency of the detec- 
tors used was insufficient to provide positive results. The 
lattice transfer function was electronically simulated and 
analyzed by the system confirming the reliability of the 
equipment. E 


Several methods are suggested for increasing detector 
efficiencies to provide signal levels which will permit 
extraction of information about the kinetic behavior of the 
test lattice. 
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CHAPTER 1 
INTRODUCTION 


l.l Background 


In studying the neutron characteristics of any nuclear 
peactor core assembly it is desirable that the research be 
done a8 economically as possible, and preferably that 
several independent methods be used to determine the 
desired information. “he MIT Lattice Researoh Project has 
for some time Studied ,arious reactor parameters in a heavy 
water-moderated exponential facility located adjacent to 
the MIT Reactor (1-5) See Figure l.l for a cut-away View 
of the MIT Lattice Research Facility. 

In recent years, considerable use has been made of 
the pulsed neutron technique in determining kos for test 
lattices in order to provide data supplementing the more 
common static parameter ری ”سے ہک‎ 097 It is the purpose 
of this thesis to explore the utility of the noise analysis 
technique as a means of complementing the pulsed neutron 
and other methods of measuring the various parameters 


measured in the MIT heavy water lattices. 


1.2 Noise 

The term "noise" refers to fhe random or stochastic 
fluctuations in the output signal from a system. Reactor 
noise is thus the random fluctuation of power level due 
to the discrete n ture of the multiplication process. 
These fluctuations are generally small. and urually only 
about one percent of the total power bə In order to 


observe this noise, a detector must be used. and often 
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the term reactor noise is used to describe the fluctuations 





of the output signal of the detector. It is important to 
remember, however, that the detector only seeys a token 
sample of these fluctuations a.d that, owing to the finite 
size of the detector, the spatial-dependence of the noise 
measurement is blurred. Only recentiy has the space-depen- 
dence of reactor noise been closely lea e? Throughout 
this paper, however, spaoe-independence will be assumed 


since this will prove alequate fO preset PULpOe sd o 


1.3 Summary of Theoretical Curaiderations 

Since a reactor depends upon neutron multiplication for 
its operation, it is reasenable to expeet thgt among the 
various neutrons which are registered by a detector, many which 
have a common ancestry will produce pulses superimposed upon 
the d.c. level of the signal. The number of such related 
pulses will depend upon the detector affiulenty as will be 
seen later. While the detector d.c. level is itself main- 
tained by the accumulation of many pulses, we are unable to 
extract any information involving the discreteness of the 
multiplication process from it, and it thus holds no interest 
D aus here. Thus one can ineteed focus attention upon a 
fluctuation in the signal due to the daughter neutron of a 
given fission, followed a certain 2: later by another 
impulse due to a grandaughter neutron of the same fission. 
Since this time interval At is a fuüciilon of thé neutron 
lifetime in a given system, it sod be ına)? to extract 


information about the multiplyine eres s of thə medium 








a 
by submitting the ac (òr abide) portion of the detertor 
signal to analysis. One method of analysis is to determine 
the autocorrelation funétion, which is defined as the expece- 
tation value of the Lrpduct nf tlc signal at a Lime t, and 


the signal at a time t» s 


ÉD) = Efalt paita) a (1.1) 
where for simplicifv fc signal if Walen Го be the neutron 
density n(t), and 7 2 ۔؟۔ ہ+‎ Pollowing tke development of 
Bierman, et "SS ола Can @WAlnate B ut 6) hy first cone 


sidering the standard reastor kinetits equations in their 


spaceeindependent, linear, timsB- invariant forms 


dn — CMS (1.2) 


dt 3 
dc B. 
Bı جا‎ | 
"at = À 4° i + K 2 n 9 ( de o 3) 


Since only premit n*utroón hema is DI importance; 


these equations may he put in the forms 


dn لمت‎ 08 = ۱ | 
q£ + on "udis ng * 8 T(E} (1.4) 


where @ -—... the promp^ neutron decay constant, and 
the inhomogeneous terms are бсоселсогеа ав Фе forcing func- 
ölən f(t). 

The solution to Eq. (1.4) may be written a8 a convolu- 


tion integrals 


n(t) = A ۳ Los 
2? 900 : 





5 
where W(t-f) is the impulse response function of the system, 


This function i8 obtalned by normalizing the homogeneous 
Solution to Eq. (1.4) and i8 found to bes 


a exp[-d(t-t]] t> 


W( t=T) = (1 6) 
645 
One may evaluate Eq. (1 1) giving: 
a =a? 
Pap ET) = 2 e o (1 уд 


If one were to assume that the forcing function were 
instead a burst of neutrons, aB in a pulsed neutron ехрегјә 
ment, f(t) could be approximated by a Dirac delta ó(t) giving 


as a solution to Eq. (1.4): 

n(t) = ae OF 

At this point, the complimentary nature of the pulsed 
neutron and noise analysis experiments becomes evident, for 
each provides an independent method for determining the 
prompt neutron decay constant a. 

Various noise analysis methods are available for deter-= 
mining a by studying the a.c. detector signal in the time 
dekain, 52? however, by transforming Eq. (1.5) into the fre- 
quency domain by the use of Fourier transforms, a method of 
analysis which involves less data processing can be employed. 

The Fourier transform of the autocorrelation function 
is called the power=spectral-density function and is given 
bys 


1 eo | e At | 8 
Hoa) =F] un(e ۱ (1.8) 


= 





which when combined with Es. (1 Y). Liveza 


ra 
۳ = m 
Wan ا کے‎ | 1109) 


57.2 

The spectral shige af Djs bi ај ie iho: AE p Lew pass 
filter with a breakpoint at a = '. Immun, Y: a Fourier analy- 
sis of the detector eignal Ya perforesa, if Should be possi~ 


le to determine oy the promgr Arbre dena) «“onatarto 


1.4 Detector Eft! jc, 

In oräücı bu com. Ne outil 1 fr ت2ا ۴چ ت۹0‎ 
and noise measurement, “ism EY ai ' 1 neped that both 
experiments are hampered С̧у other "imiricus signal levels. In 
the pulsed neutron expert bien hese bare the form of residual 
neutron levels due to Һи neutron ۴ in the System, and 
can be overcome by in^veesing ۴ tensity. In noise 
experiments, as we Shall дә An Unie seorion, the detection 
process itself provides a mall” hacuwr3vpBd power spectrum, 
which can only be overoume by длр 2 лн detector efficiency 
or by performing сгрбә-ооутР 1061 Һи әсл (тб. These limita- 
tions can be deseriped Oy ۳۵۱ Tle oáanberound as an 


additive constant tê burn 80۳٣۷۷٣۷۷۷٣ 








E 
t) = ane 1.10 
N(t) Ар + B. ( ) 
Ba ) 


where now a figure of merit {nr ۳ system can be defined ass 








T 
utility and feasibility nf noise experiments thus depends 







upon the magnitude of the figure of merit 

By pursuing the analvtical method of Cohn 9) , and extend- 
ing it to the case of a subcritical reactor. the effects of 
"detector efficiency can be studied. The basis of this method 
rests on what 15 known as the Sohottky formula, which was 
originally developed tw jescrihe the noise in a temperature 


limited electronic 2122 IL is: 
бу 
Jbu - 7 , (1.12) 


where qi» is the spectral density cf mhe diode current 
noise in amps” sec, e is the charge sarried by each electron 


in coulombs, and m is the average number of electrons flowing 


per second. Cohn adapta this Pormula to the reactor case 
giving: 
7. “4 
ام4۱5‎ 3-2 | qi m, (1.13) 
B 


where «18, I "S is the spectral density of the noise equiva- 


lent source in neut nans? sec 5, A, te the net number of neu- 
trons produced in a nuclear process of type i, and m, is the 
average number of nuclear processes of type i occurring per 
second in the reactor. As shown in Appendix A thís formula- 
tion results in the following expression for the noise equiva- 


lent source in a subcritical agsemoLy: 


^ 


- c 1 
ав.129 = 22 E 22 Е 4 Еј? 2р4а1)) (44) 
0 — 7 1 i 
This source is transformed ^v ” c£ reajfor 2”. > transfer 


— 


function T(w) into the neutron speoir sity: 
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cinçə)1”? - i7(o) Pagi (1.15) 


By considering the power sredtral density of the detec- 





tor current due to both the correlatud aut uncorrelated neu- 
trons, it can be shown bka tike Ua Aithor current auto= 


poner-spectral-density 16: 


"ə 


2 ok. 1 ir ә 
«nr? - 2Q SET ҝ —-= UE» —. Sri mn -2Y +1) o , 
ҝ| T y y GO = 
(1.16) 
where 
fissions detected 
= 2 Ké m. + Z ee رس‎ 
€ detector efficiency FTastons occurring 
kon 
Бә = fissiun mle = -= 
u^ 
k  k(1-6) 
о = (1-8) 
yp = average number nf neutrous emitted per fission 
О = charge collected per fienien detected 
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"XE 


Thus, for an autotorrela"lon meaaursmént, the detector 


-3 
gə 
c 
< 
To 
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efficiency must be such that the rigbt-hand term is greater 
than unity, or for c << Gà 


۳ اد‎ 
gö U A 


k k 
P : Bons 2 7 + 1) 


y) 


> ٦ 





> 


Using the accepted RD of 2.5 for the factor y and 


7.47 for ve, this bet 05: 


(1.18) 











Figure 1.2 displays thie mesnmiorien for k "о 90. 

Since detector effibiíiencieus above 1072 are very difficult 
to obtain without significantly ularurbiue the system under 
study, a major stumbling hime To be wen of noise analysis 
in far-subcritisal svsteme {IN DO} hevomes obvious. Some of 
this problem could be alleyteted if the uncorrelated noise is 
filtered out. By ual e бәс Anhretanrs aml multiplying the a.c. 
components of their ser. rheeg pulse wRnich did not coin- 
cide would tend to ^$nec? per meetéeds fron which did coincide 
would reinforce eacb olné- smi, own 4 vir ef time, the 
correlated portion cf the spectrum would tend to be emphasized. 
This method of cross-correjlation wac used by Seifritz, 
Stegemann, and Vatb in making wea3urceaen*ts on the coupled fast- 
thermal reactor STARK LI) Pheip césultf& indicated that with 
their equipment, a ratio cf рәфгбллеЗ to uncorrelated noise 
of O.l was sufficient to giye meApiwgüu! resulta if croas- 


correlation techniques were uial., This ۱8 a factor of twenty 


improvement over the ratio of z «mich Lt) needed in auto 
correlation experiments. 07 course, tbe relatíve advantage of 


cross=correlation over aut -correlation experiments depends 


upon the individual amaly®! \ulfew And how the signal is 
analyzed. Figure 1 1 3howr I! ""us-gervrelated detector 
efficiencies required based 00 +32 ۴ - Notice that 
detector efficienties gf tivese br nne pepcent are still 
required fora k of JU | 3, whieh 13 the range of pri- 


mary interest in MIT ۔‎ ۴ n». 





This thesia if none rue | jx: r Aanembliink a 


reliable analysis system to penptuew! ۵ auto* and 
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11 
cross-correlation speoTy/üal imyo surpudubu. From the start, it 
was realized that ond such ¥ swat were construehed, the 


Problem of attaininw suftielsuti, larg üstector efficiencies 





would have te be anived. 

Detector efficiency ig ۷۱ a matter of detector size 
and neutron sensitivity ОЈ all déleotope readily available, 
the Westinghouse Tute iT ionhestion cwamber provided the 
best combination cl neuzur deanedlayity (4 x ic^ amperes/ 
neutron/cm”/sec ) ang siz, (Va pw. ct şgansıbive area). -As 
such, this detertur war 4311. oa ş$finienk ensuğgh (in the 
hohlraum and lattice, * = |) s 10 |} Boose asure the spectral 
density; however it served tee murroges Fivat. 16 provided 
a noise signal corparalle iv side fo those whach would be 
studied in the lattice latev, Jha thus provided us with infor- 
mation on tne design ral rara For the analysis system. 
Second, any very large ıesgusnsee In the noahlraum which might 
be present due tao spatial wlfechs men still be detected. 
Thus, both auto- and i roesascorrm 400.9 menmmurzments were made 
in the lattice to insure that a reliable system would be 
available in the future. concurrently, a separate study was 


being made of methods fcr in.) lms d&tector efficiencies 


to the levels needed, (12/ 


Chapter 2 describes ine wort dite che fiiitial phase 
3f Ey^AME sf3embly ard Yugdi2^us 0 sura. Chapters: 3 
ATI A i : اھ‎ - by both 
le zd - hapter 5 





12 
CHAPTER 2 


ASSEMBLY OF EQUIPMENT AND PRELIMINARY 
HOHLRAUM MEASUREMENTS 


2.1 Assembly of Equipment 


Several methods of analyzing noise signals have been 
developed, often classified as to whether the analysis is 
done concurrently with reactor operation, or done afterward 
by analyzing either an electronically or mechanically (pen- 
recorder) recorded dignas. 2) Whether the analysis is done 
digitally or continuously, further specifies the method of 
analysis to be used. The choice of method of course depends 
upon financial and/or time limitations, and also upon the 
degree of accuracy desired. 

For present purposes, it was decided that the use of 
analog computing equipment would provide the most flexible 
system of analysis. Not only would it provide a real-time 
method of continuous analysis but it would also be adaptable 
to other uses such as reactor dynamics simulation. For this 
reason one George A. Philbrick MU/DV Duplex Multiplier/Divider 
and three G.A.P. Model HK Operational Manifolds (each con- 
taining five chopper-stabilized operational amplifiers), were 
acquired. The latter afford the capability of performing a 
total of fifteen separate linear operations on one or more 
signals, whereas the former performs the non-linear operations 
of multiplying or dividing. (See Appendix D for a brief 
discussion of the characteristics of operational amplifiers 
and a description of the basic ۵ ` aged in performing 


the linear operations.) The method of a + jsis used was to 








43 


take the Fourier integral of the sav e-cermelae mm function 
with zero delay time. A block diagram uf the cimmuttry is 
shown in Figure 2.1. 

After using the heterodyntog tschnigqu- cf multiplying 
the detector noise by a Bine wave of known frequency, a low 
frequency band-pass filter or fixed barovıith A- is used to 
pass only those nolüe compeoneptbaü wnoaé feequeney differs from 


that of the oscillator hv Jean than the Yiitoo bandvidth, 


The resulting sianaı is aw 00 1 1 ۰۷ 1 1 «3 by 
the integration time ivan ۷ 0 Ел 'pectral 


density at a frequency co, - A, in a bandwidth AD. 
C^ 


Scanning the spectrim is au. iipiiwhed by vary lug the oscilla- 


= 


tor frequency w , over the ines tg 


2.2 Preliminary Experimentat in 


Initial assembly of tl -quipment od niit hiton 
were done by the author tos 7! (th | | = 113) 


As stated earlier, this first pha®w whs te serte WO purposes. 
First to determine if the aniy fine 1 | аз performing 
correctly, and secondly to determine if the rhermal column 


and hohlraum served as a trüe white noise source of neutrons, 


This latter question is inp " singe, Wa ofm be seen by 
Eq. (1.15); determination of the lattide banat er function 
requires foreknowledge of the m Tere i P the source, 


Only when the external source is 1tsclr pandom may the 
2 cx Ñ 
(110) 


simplified pile noise formula fui 1 
For these measurements, í dh صن‎ iel 6377 loniza- 
tion chamber, was placei in 11 ah a hole bored 


through the graphite wall in: | 
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13 
Preliminary measurements showed that the total d.c. 


detector output was 175 pamps with a noise component of 
about 1 namp. Since it was desirable to have a voltage 
input to the vacuum-tube equipment on the order of one-half 
volt or larger, the need for a high gain preamplifier became 
apparent. To avoid excessive 60 cycle pickup and to intro- 
duce as little instrument noise as possible, it was decided 
to use an Analog Device Model 106 Solid State А.С. орега= 
tional amplifier, powered by four F4BP 6 volt batteries. In 
addition to meeting technical requirements, the amplifier 
has the added advantage of economy costing $26, For similar 
reasons, the detector high voltage was supplied by 10 67.5 V 
batteries. To avoid saturating the preamplifier with the 
d.c. level of the incoming signal while still obtaining maxi- 
mum amplification of the noise component, a buckling voltage 
E by a potentiometer was inserted just before the 
input. Also, since there was a large low frequency variation 
in the neutron flux, apparently due to control rod motion In 
the reactor, a high pass network (R, and C.) with a break- 
point at 2.5 cps was introduced (see Figure 2.2). Thus the 
preamplifier response is: | 
+ ee), 
$ Sos an 
(2R, * R,)(1 + PR. T R.) 


(2.1) 


which for large s = j¥ becomes 


SM 


— Á—ó[ R Ы 
Lan 1 


Although a high frequency roll-off due to finite space charge 
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R, = 510 KQ 
Co = 220 uf 
R, = 100 KA 
P, -500 KR 
E = -12 V 
A = 150,000 
lj ————— — ہے‎ z = کے مہ ہے‎ 
E > Rz 
PUTAT. 
E G 
2 
Pi = 


Fig. 2.2 High Gain A.C. Preamplifier Circuit 
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collection times was observed, the preamp-detector frequency 
response was constant up to 5,000 cps for a detector voltage 
of 720 V. 

The noise amplitude at the preamp output was about one- 
half volt, which was sufficiently large compared to the 
0,01 V self-noise of the preamp. Two stage amplification with 
a total gain of 200 provided a 100 V peak-to-peak signal at 
the input to the heterodyning multiplier. Here, the signal 
is multiplied by a 100 V peak-to-peak sine wave of known 
frequency. Since the output of the multiplier is scaled by 
100 V, the resultant signal is 100 V p-p. The combination of 
the two filters produced a 0.88 cps pass-band centered at 
0.128 eps, The inclusion of an amplifier with a gain of 20 
between filters maintained the signal at a sufficiently high 
level to make the effects of component drift at this stage 
negligible. After squaring, the filtered signal was integra- 
ted using an averaging or "leaky" integrator, with a time 
constant of 100 seconds. The circuit at this point is shown 
in Figure 2.3. Results of the preliminary measurements in 
the hohlraum are presented in Figure 2.4. Two features are 
important. The first, already described, is the high 
frequency roll-off as a function of detector voltage. The 
second is the low frequency erratic peaking in amplitude. 

To investigate the latter, the effect of the reactor 
Source neutrons was considered. Since these neutrons orig- 
inate in the critical MIT Reactor, a low frequency breakpoint 
in the reactor neutron spectral derg?” would be expected at 


w=p/L . This value is approximately .0075/.0015 for the 
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MIT Reactor, resulting in a break at f =< /2r = 0.796 cps. 

In order to determine the amount of correlation intro- 
duced into the lattice by these source neutrons, one can 
assume that the lattice is a detector and apply Eq. (1.16). 
To do this, the efficiency of the lattice for reactor کے‎ 
trons must be estimated. Less than 0.02% of eno cea 
neutrons pass the hohlraum. (44) This value may be taken as 
an upper limit for the efficiency of the lattice. Using 
Eq. (1.16) for the case of a critical reactor (k= 1), with 
a = 4,8 sec 1, it is seen that above about 5 cps, correlation 
of neutrons in the lattice should be negligible. Since the 
frequency range of interest in the lattice is about 3,000 
cycles, the effect of the low frequency correlation may be 
ignored. Thus the hohlraum may be considered to be a white 
noise sourc# of neutrons for the lattice. 

For this reason, ée behavior of the spectral density 


curve of Figure 2.4 below 100 cps caused some concern. 


2.3 First Phase Conclusions 

In order to check on the reliability of the system, 
measurements of the spectral density were taken several times 
on different days. Above 50 cps, the results were found to 
be reproduceable within the predicted experimental error (see 
Appendix Di. Scattering of data points below 50 cps cast 
doubt upon the reliability of the system at low frequencies. 

Preliminary investigations were concluded at this point 
with the next objectives being to? 


a) attempt to improve detector efr.: iency 
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b) investigate more thoroughly the lov frequeney 
. response of the system 


c) investigate frequency response in the lattice 
itself. 


The problems of developing a detector with increased effici- 
ency for lattice noise analysis is the subject of another SM 
thesis 1?) by M.G. Johnson, currently underway. The next 
chapter describes work done by the present author on the 


items (b) and (c). 





ee 
CHAPTER 3 


LOW FREQUENCY TESTS AND AUTOCORRELATION MEASUREMENTS 


Section 2.3 deseribed the erratic behavior of the analyzer 
at low frequencies. It was necessary to determine the cause of 
this behavior before proceeding to lattice autocorrelation 


measurements. 


3.1 Introduction 

Up to this point, the oscillator signal had been provided 
by first a General Radio Unit R-C Oscillator Type 1210-C and 
later a Hevlett Packard Test Oscillator, Model 650A, the latter 
providing the signals of lowest frequency down to 10 cps. In 
order to look more closely at the low frequency response, a 
Waveforms Model 403B oscillator was used, providing signals as 
lov as 1 cps. The filters were also modified as shown in 
Figure 3.1. 

This time the spectral density of the hohlraum noise was 
investigated in the range 1 to 200 cps, with the results shown 
in Figure 3.2. As can be seen, results were again erratic, 
with a general trend of higher amplitudes at lower frequencies. 

The erratic low frequency response suggested the possibil- 
ity of component drift, while the increase in amplitude at 
low frequencies suggested that the d.c. portion of the detec- 
tor signal was not being completely cancelled, allowing the 
carrier wave to be passed to the filters. Referring to the 
latter problem, note that a one-half volt d.c. signal at the 
noise input to the multiplier, when multiplied by a 100 Y 


peak-to-peak carrier signal, of frequen." ~ |, will result in 
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a one-half volt peak-to-peak sine wave at the multiplier 
output. At high frequencies (such that u , 18 much greater 
than the center frequency of the filter) this component is 
filtered out leaving only the signal due to the a.c. (noise) 
pertion of the signal. However, since the filters are not 
square (effectively 24 db/octave attenuation as seen at the 
averaging circuit input), at low frequencies this sine wave 
would not be attenuated enough to be negligible. Total cancel- 
lation of the d.c. level of the Signal was insured by adding 
a blocking capacitor in series with the first amplifier. 

Multiplier drift (rated at 419) which was measured to 
be as much as 2 V in one hour, accounted for the erratic 
behavior at low frequencies, This drift was minimized by 
adding a ventilating fan to provide better temperature stabil- 
ity; however, drifts of up to 0.75 V per hour were still 
encountered. Figure 3.3 shows how these measures served to 


reduce these problems. 


3.2 Simulation of Lattice Measurement 

To insure that the equipment was indeed operating 
properly, it was decided to simulate a lattice measurement 
by adding a capacitor to the preamplifier, in parallel with 
the feedback resistor. Neglecting the low frequency filter 


network, the preamplifier has a gains 


e R 
مھ"‎ 9 0 (3.1) 
SUITES. 


With R, = 510K (1$) and C, 2 .20. nf, ^"sagured) this results 


f 
in a simulated breakpoint at: 
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er 
1 
f = ee = 314 + 4 eps 
2 HAC. 


At the same time, the filters were modified to widen the 
pass band, which would allow shorter integration times (see 
Appendix D). The filter modifications are shown in Fig. 3.4, 

Results of the Simulated lattice measurement are shown 
in Fig. 3.5. The breakpoint occurs at that frequency corres- 
ponding to a 6 db decrease in the squared signal and is 
estimated to be 320 + 20 cps, which is in good agreement with 
the predicted value of 514 b X (5, 

3.3 Lattice Autocorrelation Measurements 

Having been satisfied that the system was now functioning 
properly. attention was turned to making a test run in a 
D,O-moderated, uranium-fueled lattice. Estimates of Kor? and 
6۵ for the lattice installed in the exponential] facility during 
this thesis work are developed in Appendix C, indicating that 
a detector efficiency of around 302 would be required to 
enable determination of a for this particular lattice.  None- 
theless, since it would be Several months at least until a 
more reactive lattice was scheduled for installation, it was 
decided to make a test run to insur. that the equipment would 
be ready for use, These measurements were made by Suspending 
the 6377 boronlined detector in a lead-weighted aluminum tube 
in the center of the lattice core, In order to do this, the 
Seven center fuel rods had to be removed from the 211 rod 
lattice. Since the multiplication facter of the assembly was 
already 30 low as to preclude measur '^ of kinetic effects, 


the small additional decrease in reactiv =. caused by removing 
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the rods was not considered significant, First attempts at 


taking measurements revealed that the 75 foot cable from the 
detector to the preamp caused excessive oscillation in the 
preamplifier. Although this oscillation could be eliminated 
by use of an RC compensating circuit, the signal strength was 
then so small as to be comparable to the preamp self-noise. 
Therefore, the preamplifier was placed on the lattice top, 
eliminating the need for both the long input cables and the 
compensating network, Such an arrangement produced lattice 
noise signals five times greater than the preamp self noise. 
With the equipment set up in this manner, the results shown 
in Figure 3.6 were obtained. Again, the white-noise spectrum 
due to the low detector efficiency is evident. 

Мо further work was performed on autocorrelation measure- 
ments because the cross-correlation technique to be described 


in the next chapter appeared to hold more promise, 
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CHAPTER 4 


CROSS=CORRELATION MEASUREMENTS 


pl] ۳ 


As mentioned in Chapter 1, crossecorrelating the 
signals should decrease the detector efficiencies needed 
for noise measurements, The cross-correlation function 
is defined as the expectation value of two independent 


signals: 


E İg (6 əsə) (4.1) 
Although 2, and H. consist of noise signals super- 


imposed on d.c. levels, we need consider only the noises: 
Ø lt) =u (t) + c(t) (4.2) 
gəc) = U,(t} + Colt) 9 (4.3) 
where U and C are uncorrelated and correlated noise, 
respectively. 


The cross-correlation function of these two signals 
is given bys 
EİZ, (6322(6)İ - 
Efu (6)02(6)) + E[0,(6)02(5)] 4 ејб (ес (6) || 4 Е|2, (2249) 
(4.4) 
Since the expectation value of the multiple of two 


uncorrelated functions is zero, all uncorrelated terms 


disappear and: 


2 ]2 (OyətƏl - si ztel < (4.5) 


Thus, corresponding to Eq. 1.16 for auto-correlation 


experiments, one can write for cross-correlation experiments: 
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E (4.6) 


2 ә ۲ 
۱1081۱ eg 
cə + G 





where now € , and Ед refer to the efficiencies of the 
separate detectors. This would seem to imply that measure» 
ments could be made independent of detector efficiency, as 
long as the analyzer output signal is sufficiently large. 
However, as stated earlier, Seiffritz, et al. found that 
cross-correlation decreased their efficiency needs only by a 
factor of 20. This is due to the fact that the cross-terms 
in Eq. (4.4) are not totally uncorrelated. A factor of 20 
reduction in efficiency requirements, however, is a highly 
desirable, if not essential, improvement if one is to make 


noise measurements in subcritical systems. 


4.2 Equipment 


The cross-correlation measurements required tbe addit- 
ional equipment shown within the dotted lines in Fig. 4.1, 
consisting of one additional multiplier (G.A.P. + 25 V MU/DV 
Duplex Multiplier/Divider), 6 operational amplifiers, 
another preamplifier, and an additional Westinghouse Type 
6377 boron=lined ionization chamber. 

The two detectors were suspended in the high-flux 
region of the lattice described earlier, side-by-side, in 
Separate lead-weighted aluminum tubes. Again the preampli- 
fiers were placed on the lattice top, as close to the detec- 
tors as possible to avoid preamp oscillations. Both detec- 


tors were operated at 670 V. 
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4.3 Results 
Recall that for autocorrelation measurements the filtered 
signal is squared and then averaged. The squaring process 
naturally results in a large hə ə: whose average 
will also then be non=2 ie orossecorrelation measurement, 
in that it eliminates the uncorrelated noise signals, should 
produce a signal whose magnitude relative to the auto-correla- 
tion signal is given by the ratio of Eq. (4.6) to Eq. (1.16). 
For very low efficiencies, this ratio is merely the right-hand 
term of Eq. (1.16) and iss 
“al” 


(II, |I» 


1 Ke. (4.6) 





For €x10"/, and for the lattice under study: 


2 


-6 


cit, | 
^ 5 SNO x VU 


«IIl 


Thus since the signal level in the auto-correlation measure= 


256 100 V, the cross-correlation output signal 


پ - 


ments was Cl L.I 
must be «1,17» SS e This signal is so low as to be 
overridden by equipment drift. While observing this drop in 
net output confirms predictions as to the effect of converting 
from auto- to cross-correlation, the low signal made it 
impossible to extract a spectrum. Even with crossecorrelation 
techniques, therefore, higher lattice multiplications and 
detector efficiencies will be required. 


Experimentation was concluded at this point. Conclu- 


sions and recommendations follow in Ch er 5: 
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CHAPTER 5 


SUMMARY, CONCLUSIONS, AND RECOMMENDATIONS 


k. $t ear 


5.1 Summary aha’ Conclusions 


A reliable system for making auto- and cross-power noise 


spectral density measurements on a subcritical nuclear 


assembly was constructed. The system was demonstrated 


capable of providing signal analysis over a frequency range 


of 1 to 5,000 cps. 


The assembled system waa tested in several applications: 


a) a signal simulating a reactor noise spectrum 


b) 


was analyzed to demonstrate successful system 
performances under optimum operating condi- 
tions (Chapter 3, Section 2). The measured 
breakpoint frequency 320 +20 sps agreed quite 
well with the known input of 314 + 4 eps. 


The MIT exponential facility hohlraum was 
studied using the auto-correlation technique. 
None of the results (Chapter 3, Section 1) gave 
reason to doubt the assumption that the hohlraum 
provides a white source of neutrons to the test 


lattices. 


Both auto= and cross-correlation measurements 
were made in a lattice. In both cases it was 
shown that lattice multiplication and detec- 
tion efficiencies must be increased before 


noise analysis can be used *^ extract neutron 


decay constants (Chapters 3 ani 4). 
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Since detector efficiency requirements vary inversely 
as the smuare of the multiplication factor of the system, 
future studies should concentrate on lattices with Koff as 
close to the present lattice safeguard limit of Koff x 0.9 
as possible. Since cross-correlation of two signals 
significantly reduces the efficiency requirements, this 


method provides the greatest promise for successful analysis 


at reasonable efficiency levels. 


5.e Recommendations 

In order to attain the required detector efficiencies, 
the utility of gamma detectors in the lattice annulus 
should be fully explored. Such detectors wouid not disturb 
the multiplying system, but would detect both prompt fission 
product gammas, which escaped from the lattice core, and 
those gammas produced by neutron capture in the cadmium 
sheath surrounding the lattice. Alternatively it may be 
possible to develop a detection system which detects the 
Cerenkov radiation produced in the moderating D50. Further, 
the 2% enriched UO, lattices, scheduled for study beginning 
in February 1967, should provide the capability of attaln- 
ing Karf 7 9:95 which would also reduce detector efficiency 
requirements considerably. 

Finally, in future work, where large volume detectors 
must be used, the space dependence of reactor noise measure- 


ments should be carefully considered. 
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APPENDIX A 
NOISE-EQUIVALENT SOURCE FOR A SUBCRITICAL ASSEMBLY 


Neutron processes in a multiplying system are analogous 
to the random flow of electrons in a diode. Using this 
Similarity, Cohn has adapted the diode noise equation of 
Schottky to the noise in a critical reactor. 10) my simple 
modifications it is possible to extend this treatment to the 
case of a subcritical system. 


Cobn's adaptation of the Schottky formula is: 


eS E 2 2 q, f, 9 (A.1) 


<1s 
z i 


where: «i18, 17» is the noise-equivalent source; 
q; = net number of neutrons produced in the nuclear 
process of type i; and 


average number of processes of type i occurring 


81 
hr 
il 


per second in the system. 
Table A.l lists the processes, their rate of occurrence, 


and the net number of neutrons produced in that process. 


TABLE A.l 
CONTRIBUTIONS TO PILE NOISE SOURCE 


Average Rate Net Number of 


Nature of Process of Occurrence Neutrons Produced 











Nonproductive leak- 57 : e SR 
age or absorption P. 

Fission giving rise n F " 
to N prompt neutrons AA + F EN 02 


External source provid- 3 
ing neutrons to main- — (1- k.) 1 
tain the pile 4 : 


LI‏ سے 
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In this formulation delayed neutren effects have been neglec- 
Bed. Here 


A macroscopic cross-section for ali non-productive 


processes including leakage, 


P 


macroscopic cross-section for fission, 


n = total number of neutrons in the system, 


A =/k = neutron generation time in sec. 


prompt neutron lifetime in sec, and 





A and F are subject to the neutron balance soudition 


; 


d 


۱ 2 5 
a Ҹ̧ (A.2) 





where » is the average number of neutruns produced per 

fission. P. is the probability that N prompt neutrons will 

be produced in any one fission, subject to the conditions 
00 


En = 1 (A.3)‏ م2 
N=1‏ 


Co исы. SA 
ZN Py = (1-8) Y = 
NEST 


(A.B) 


4 


00 
2 
AN Eu 


y? 


3 
A 


- (1-8)? v^ « ; (A.5) 


N= 
Substitution into Equation (A1.1) yields 


00 
IE. = JF] A+kF “ (N-1) Py + (1 ul ə (A.6) 
Applying conditions (A.3), (A.4) and (A.5) yields 


i 28 
وگ واه‎ «22) |3 - E| 4| £-t*-22 41)| + [1=k]0. 
o 7 زر‎ ۱ 


(A.T) 


The first term in square brackets is tr: term due to 





ho 
nonproductive (e.g., (n,y) capture and leakage) processes; 
the second is the term due to ffssion processes, the third 
is the term due to the external source. Note that the exter- 
nal source term has been assumed to be frequency independent, 
and that the entire noise equivalent source is thus a white 


noise source. Since the fission rate ise 


Ру = =F : (4.8) 
Eq. (A.7) becomes 
237 Ko ke vm uc 
> | ام5‎ >= 9, (İz - > - ky + — - 2) 4 7 (A.9) 
The source neutrons are then correlated by the multiplying 
process. The correlated component of the noise spectrum is 


obtained by multiplying equation (A.9) by the nuclear system 


transfer functions: 


«In(2) ^». - 17(2)I* «ts, it» 


2 
Zi, k k 
1 
= —S|2 - -È = k. + By? -2y +y )=— 
K, y y WE + a 


(A Aen 
where a - Ba the prompt neutron decay constant. 
The average current I in a detector due to all fissions 
detected will bes 
m 


= €Q => = EQF, (A.11) 


where Q is the charge collected per fission detected. Super- 
imposed on this d.c. level are noise components due to the 
correlated neutrons and the randomne- ОТ the collection 


process. The former is given by Eqs. (A.10) and (A.11) 
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2 2 ax ^ 2 
€I1 |^ 52€ ۱ء مقر‎ > (A.12 
P O) ) 
DER k.. گا‎ 
=—— s |2 - Ë - k. + ——(v“ - 25 ges. 
07 У ۹-0 i 
(A.13) 


The latter can be obtained by using the Schottky formula: 


«4x,1^» - aq”er, . (A.14) 


For € € € 1, these two terms are first order uncorrelated 


and are added in quadrature to give 


€x k к 2 
2 o 2 _- l 
LIT (w) ("> = 2Q°CR 1+ 2چ‎ - = - k. + y^-2y +1 əə 
a? al = y p T er 


(A.15) 
This result enables estimation of detector efficiencies 
needed to extract information from noise signals as shown 
in Chapter 1, Section 4. 

Note that throughout, it has been assumed that the 
analysis equipment did not introduce any spurious frequency 
dependence. Such an assumption is justified by the “white" 
noise spectrum displayed in Figure 3.3. 
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APPENDIX B 


OPERATIONAL AMPLIFIERS AND ANALOG CIRCUITS 


An operational amplifier is a high gain amplifier with 
certain input and output impedances as shown in Figure BL. 


Analysis of the circuit shovs that: 





e 2 
o Б 1 
Ha ۱ (хәз 
Dg 1) 


where Zr feedback impedance 


2 input impedance 


i 
A = open loop gain. 


7 
Note that if A > O (ZE 4.1), this reduces to: 
i 


4 
E 28 (B.2) 


(D 


P 
e 


Since Ze and Z; are complex quantities, various combinations 
Of resistances and capacitances may be used to perform mathe- 
matical operations. Table B.l illustrates the circuits 
which are used to perform the linear operations in the 


analyzing equipment. 





Fig. B.1 Operational Amplifier Cr rcuit 
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APPENDIX C 
ESTIMATION OF Keff AND o 


Estimates of Kas and a can be made by using the 
THERMØS computer code .(1*) his integral-transport theory 
code, based on the Wigner-Seitz unit cell model, permits 
the calculation 0۶ و‎ 
the thermal utilization 


7 the average number of neutrons produced 
per thermal neutron absorbed in the fuel 


22 the average macroscopic absorption cross- 
section in the cell 


Combining these values with reasonable estimates 


ofs 
T = age to thermal 
& = fast fission factor 
p = resonance escape probability 
D = diffusion constant 
Ba = geometric buckling, 
Koo can be estimated by: 
ü (n€ pf )exp(-B,” 0) | nai 
1 + LB,” 
where L = = (diffusion length)“ : 


o "ül 


The lattice investigated in the present studies had 
the following specificationss$ | 
Pitch 2.25 in. 
Enrichment 0.947% y=? 


Fuel Density 18.90 gm/ce 
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D50 Purity 99 .44 MPC 


The unit cell dimensions (outer radius minus inner radius) 


are: 

Puel 0.49149 cm 
Air gap O .00762 cm 
Cladding 0.08636 cm 
Moderator 2.41491 cm 


Using this data, and the procedures outlined in reference 


(14), the following estimates were derived from the THERMAS 


output: 
3, = 0.007583 cm“! 
f = 0.969 
e = 1.447 


For this lattice, good estimates for the other quanti- 
ties (9) ares 


D = 8.75 cm 

p = 0.92 

C < 1.01 

? * 115.4 em” 
B, = 0.0031 cm”? 


These values result in: 


~d 
ka, 1.30 


Koff = 0.65 


£ 


8 


4,35 x 1072 sec 


Be 


l = k 
o — ƏT -v805 sec”) 


L 
Consulting Pigure 1.2 one sees itn. , in this lattice, 


required efficiencies for successful noise analysis experi- 





ments vould be: 
for autocorrelation: 


for cross-correlation 


30% 


1.5% . 
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APPENDIX D 


COMMENTS ON OPERATION OF EQUIPMENT 
р.1 Equipment 
The equipment assembled to perform the auto spectral 
analysis consists of the following (see Figs. D.1 and D.2). 


1 George A. Philbrick (GAP) MU/DV Duplex Multiplier — 
Divider - + 100 V input-output, Ser. 0297 


GAP Model HK Operational Manifolds, Ser. 069 and 068 
GAP K2-XA Operational Amplifiers 

GAP K2-P Stabilizing Amplifiers 

GAP 2-Unit K-1 Modular Assembly Units 


~ 0 0 0 N 


Battery Powered A.C. Pre-amplifier utilizing Analog 
Devices Model 106 Operational Amplifier 


2 Associated Specialties Model 2 Power Supply, Ser. 
748 and 749 


1 Westinghouse Model 6377 Ionization Chamber 
10 No. 467 67:5 Volt Batteries 


1 Waveforms Model 403B Sine-Square Vave Generator, 
Ser. 1997 


1 Esterline Angus Graphic D.C. Ammeter Model AW, 
Ser. 85339 


Each linear operation is performed using one K2-Xa and 
one K2-P, vhich are plugged side-by-side into the top of the 
manifold. Each manifold can accommodate 5 pairs of these 
amplifiers. The K-1 Modular Assemblies are labeled to desig- 
nate which operation is performed. These assemblies are 
plugged into the face of tbe HK manifold so as to complete 
the circuitry for each choppere-stabilized operational ampli- 


fier. Each K-1 Assembly is labeled as to input, output, and 
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Fig. D.2 Photograph Showing Filter Module (Left) - 
Preamplifier (Right ) 
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ground connections. The eight operations performed are, in 


order (see Fig. 4.1): 
A. Gain of 20 amplification with D.C. blocked 
B. Gain of 10 or 20 amplification for noise signal 
C. Variable gain (up to 20) amplification for oscillator 
signal 
D. Band pass filter which may be switched to provide either 
narrow band (-0.3 eps) or wide band (~10 cps) operation 
E. Gain of Xl, X5, X10, X50, or X100 switchable amplifier 
F. Second filter similar to D 
G. Averaging circuit 
The time constant of the averaging circuit may be altered 
by using banana plugs to insert various combinations of 
parallel RC components into the jacks marked FEEDBACK. 
In order to perform the cross power spectral measurements, 
the following equipment was duplicated 
detecten 
preamplifier 
detector battery power supply 
een amplifiers and modular units to perform 
operations A, B, D, E, and F 
A 4 25 V GAP MU/DV Multiplier was borrowed from the MIT 
Mechanical Engineering Department to serve as the correlating 
multiplier. A GAP Model R-300 Power Supply was also borrowed 


to meet the increased power requirements. 


D.2 E 





In order to insure proper operatu- n, the following had 


to be observed: 
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a) Keep equipment well ventilated (e.g. by fan) 


b) Maintain largest signal possible throughout with- 
out exceeding rated signal input limitations 

с) Avoid ground loops 

d) Minimize d.c. amplification at preamp by varying 
resistances in blocking circuit 

e) Avoid preamp oscillations by keeping input cables 
to a minimum length | 

These comments supplement the individual manufacturer's 

manuals, which have been included in a booklet Noise Spectrum 


Analyzer, compiled by the author. 
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APPENDIX E 


PROBABLE ERROR ESTIMATION 


E» l Predicted Error in Amplitude 


The spectral measurements are made using an averaging cir- 
cuit with a time constant Т = R.C £ where R, and C, are the 
feedback resistance and capacitance in the averaging cir- 
cult. Using formulas developed by van der zie 645) and 
Bendat (10), the relative standard deviation in cross-power 


spectral density measurements iss 


1/2 
с” 2 + N./S + N„/S 2 (N./S)(N./S 
ے١٠ اوت‎ E 0d > 2 3 Һи )( 2/ 1 (E.1) 

5 ae T 

T 
where S = strength of signal of interest, 
Ny pus strength of undesired signals in channels 1 and 2, 
3 
Af = filter bandwidth, in cps, and 


& = averaging circuit time constant. 

If detector efficiencies are great enough to detect a 
level of correlated noise, S and N would refer to the corre= 
lated and uncorrelated noise, respectively. ۱ 

For most measurements in this work, detector efficiencies 


were so low as to permit measurement of only the uncorrelated 


signal level. For this case: 


2 
$ - şam m. "SS 


Preliminary hohlraum measurements were made with 


Af - 0.88 sec”1 


2 = 100 sec 
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Final hohlraum and lattice auto-correlation measure- 

ments were made with Af = 10 sec”? 
T = 100 sec 

These result in predicted amplitude errors of 11.9% 
and 2.5%, respectively, and are considered to be suffic- 
iently small for this work. If in the future, greater 
accuracy is considered desirable, it may be attained by 
altering the feedback impedance in the averaging circuit 


so as to increase (= گر‎ ° 


E.2 Predicted Error in Frequency 





Error in frequency measurements arise from the finite 
bandwidth of the filters and the lack of accuracy in the 
oscillator dial readings. Frequency resolution due to 
finite bandwidth is + Af, and all equipment dial readings 
are rated to within +2%. 

For example, using a filter with a bandwidth of 10 cps 
to measure the spectral density near a frequency of 300 cps 
would result in an uncertainty in frequency of 410 cps due 
to finite bandwidth and +6 cps due to dial error. This 
results in a total uncertainty in frequency of +le cps, or 
about 3.7%. This figure could be reduced by decreasing 
the filter bandwidth. AS can be seen, however, changes in 
Af affect both frequency and amplitude uncertainty, thus a 
choice of Af should be governed by desired accuracy in both. 
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